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Using the additional operators coming up at finite temperature, we calculate the 
masses and decay constants of the P wave heavy axial- vector Xbi and Xci quarkonia in 
the framework of thermal QCD sum rules. In the calculations, we take into account 
the perturbative two loop order corrections and nonperturbative effects up to 
the dimension four condensates. It is observed that the masses and decay constants 
almost remain unchanged with respect to the variation of the temperature up to 
T ~ 100 MeV, however after this point, the decay constants decrease sharply and 
approach approximately to zero at critical temperature. The decreasing in values of 
the masses is also considerable after T ~ 100 MeV. 

PACS numbers: 11.55.Hx, 14.40.Pq, ll.lO.Wx 
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I. INTRODUCTION 



Investigation of the in medium properties of heavy mesons such as bottomonium (bb) and 
charmonium (cc) are of considerable interest for hadron physics to date. These quarkonia 
play an important role in obtaining information on the restoration of the spontaneously 
broken chiral symmetry in a nuclear medium and understanding quark gluon plasma (QGP) 
as a new phase of hadronic matter. Such investigations can also provide us with substantial 
knowledge on the nonperturbative QCD and interaction of quarks and gluons with QCD 
vacuum. 

In the last twenty years, a lot of theoretical and experimental works have been devoted 
to study the behavior of the heavy mesons in medium. The J/ip suppression which can be 
considered as an indication of QGP jl| has been observed in heavy ion collisions in super 
proton synchrotron (SPS) at CERN and relativistic heavy ion collider (RHIC) at BNL. 
A plenty of theoretical works have also been dedicated to study the thermal behavior of 
hadronic parameters as well as QCD deg rees of freedom (for some of them and discussion 
on the QGP phase see for instance 



25|). 



Hadrons are formed in a region of energy very far from the perturbative region, hence to 
calculate their parameters we need to have some nonperturbative approaches. The QCD sum 
rules as one of the most attractive, applicable and powerful techniques has been in the focus 



of much attention during last 32 years. This approach at zero temperature proposed in [26 1 
and have applied to many decay channels in this period giving results in a good consistency 
with the existing experimental data as well as lattice QCD calculations. This method then 



was extended to finite temperature QCD in 27 1. There are two new aspects in this extension 



compared to the case at zero temperature 28|-|30|, namely interaction of the particles in the 



medium with the currents requiring modification of the hadronic spectral density as well 
as the breakdown of the Lorentz invariance via the choice of reference frames. Because of 
residual 0(3) symmetry at finite temperature, more operators with the same dimensions 
come out in the operator product expansion (OPE) compared to that of vacuum. 

The purpose of this paper is to calculate the masses and decay constants of the P wave 
heavy axial vector Xbi and Xci niesons in the framework of the thermal QCD sum rules. In 
our calculations, we use thermal propagator containing new non-perturbative contributions 
appearing at finite temperature, and take into account the perturbative two-loop order 
corrections to the correlation function l2a, 



31| . We use the expressions of the temperature- 



dependent energy- momentum tensor obtained via Chiral perturbation theory [32|| and lattice 
QCD js], 0] as well as temperature-dependent gluon condensates and continuum threshold 
to obtain the behavior of the masses and decay constants of these mesons in terms of 
temperature. 
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II. THERMAL QCD SUM RULE FOR P WAVE HEAVY AXIAL VECTOR 

QUARKONIA 

In order to extract the sum rules for the masses and decay constants of the heavy axial 
vector Xbi and xa mesons at finite temperature, we start considering the following two-point 
thermal correlation function: 

n^,(g,T) =zjd'x e^^-^{r{Ux)Jl{0))), (1) 

where, J^(x) =: Q{x)'y^'y5Q{x) : with Q = 6 or c is the interpolating current of heavy axial 
vector meson, T is temperature and T indicates the time ordering product. The thermal 
average of any operator O is defined as 

(0) - (2) 

where H is the QCD Hamiltonian and /3 = 1/T. 

According to the general philosophy of the QCD sum rules formalism, the above corre- 
lation function can be calculated in two different ways. Once, in terms of QCD degrees of 
freedom by the help of OPE called the theoretical or QCD side. The OPE incorporates 
the effects of the QCD vacuum through an infinite series of condensates of increasing mass 
dimensions. The second, in terms of hadronic parameters called the physical or phenomeno- 
logical side. Matching then these two representations, we find sum rules for the physical 
observables under consideration. To suppress the contribution of the higher states and con- 
tinuum, we apply Borel transformation as well as continuum subtractions. In the following, 
we calculate the correlation function in two aforesaid windows. 



A. The phenomenological side 



Technically, to obtain the physical or phenomenological side of the correlation function, 
we insert a complete set of intermediate hadronic states with the same quantum numbers 
as the interpolating current Jfj,{x) into the correlation function. After performing the four- 
integral over X and isolating the ground state contribution, we get 

where the /a and Ma are decay constant and mass of the heavy axial vector meson, respec- 
tively. The in the above equation stands for the contribution of the excited heavy axial 

vector states and continuum. In deriving the Eq. (JSj), we have defined the decay constant 
fA by the matrix element of the current between the vacuum and the mesonic state in 
the following manner: 

{Q\J,\A{q,\))=fAmAe^^\ (4) 
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where is the four-polarization vector. We have also used the summation over polarization 
vectors as 

E4'^*4'^ = -(^?"^- WO- (5) 

A 

B. The QCD side 

In QCD side, the correlation function is calculated in deep Euclidean region where -C 
~-^QCD OPE where the short or perturbative and long distance or non-perturbative 
effects are separated, i.e., 

n?.^^(g,T) = n^^,(g,T) + n;:^(g,r). (6) 

The short distance contributions are calculated using the perturbation theory, while the long 
distance contributions are expressed in terms of the thermal expectation values of the quark 
and gluon condensates as well as thermal average of the energy density coming up at finite 
temperature. 

In the rest frame of the medium for axial vector meson at rest, the correlation function 
in QCD side can be written in terms of the transverse and longitudinal components as 

n,.(g) = - 9,u^ n,(g) + ^Mq), (7) 

where the functions, n((g) and Ili{q) are found in terms of the total correlation function as 



Mq) = \q''q''n,u{q)- (8) 

Here, we would like to mention that the transverse and longitudinal components are related 
to each other, hence it is enough to use one of them to obtain the thermal sum rules for the 
physical quantities under consideration. Here, we use the function Hf(g) for this aim. It can 
be shown that this function for the fixed values of the |q|, can be written as [l^ : 

ll.{%,T)=l dq, (9) 



where Qq = — gg, and 



Pt 



(ql t) = hmlit (ql t) tanh ^, (10) 

is the spectral density. We also should stress that the function Ht(^gQ,T^ receives contribu- 
tions from both annihilation and scattering parts (for more information see |2^). However, 
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as we deal with the mesons containing quark and antiquark with the same masses, the scat- 
tering part gives zero and here we focus our attention to calculate only the annihilation 
part. 

The thermal correlation function in QCD side is obtained from Eq. ([1]) contracting out 
all quark fields via Wick's theorem. As a result, we obtain the following expression in terms 
of thermal heavy quarks propagators: 



(2vr)^ 



-Tr 



-ff,S{k)-f^S{k-q) 



(11) 



In real time thermal field theory, the function Ut (^q^ , can be expressed in 2 x 2 matrix 
representation, the elements of which depend on only one analytic function. Hence calcula- 
tion of the 11-component of such matrix is enough to get information on the dynamics of 
the corresponding two-point correlation function. The 11-component of the thermal quark 
propagator S{k) which is given as a sum of its vacuum expression and a term depending on 
the temperature is given as [33|: 

Sik) = {^k, + m)(- ^— — + 2mn{\ko\)5{e - m^)), (12) 

where n[x) = [exp{f3x) + 1]~^ is the Fermi distribution function and m is the quark mass. 
Performing the integral over ko in the q = limit, we get the imaginary part of the Utiqi, T) 
as: 



/mni(go^T) 



in" 



1 - 2n{uj) + 2n\uj) 6{qo - 2w), (13) 



where a; = y + m^. After standard calculations, we get the following expression for the 
annihilation part of the spectral density: 



Atx- 



1 - 2n 



2T^ 



(14) 



where, v{s) = \Jl — Am?/s. As we previously mentioned, we take into account also the 
perturbative two-loop order correction to the spectral density. At zero temperature, it is 



given as 



26, 



3l|: 



Pa,(s) 



37r3 



TT 
2^ 



1 + V /Tl 

2 ^2 

l + v 



+ (P^(t;) - P{v)) In -—- + Q^'iv) - Q{v) 



(15) 



where we have set v = v{s) and the functions P{v), Q{v), P^{y) and Q^iy) are given as: 
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QCD Quark perturbative g 
» = ' + — S + 
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FIG. 1. The quark propagator in the gluon background fields. 

To get the thermal version of the above two-loop order correction, we replace the strong 
coupling as by its temperature dependent lattice improved version given in 0] ( for more 
details see also j2l, Isl). 

Our final task in this section is to calculate the nonperturbative part of the thermal 
correlation function. The nonperturbative part in our case can be written in terms of 
operators up to dimension four as: 

U7{ql T) = C.l^ij) + C^^G^C^n + C,{uQu). (17) 

where, Cniq^) are as Wilson coefficients. As we also previously mentioned, at finite temper- 
ature the Lorentz invariance is broken by the choice of reference frame and new operators 
appear in the Wilson expansion above. The new four- dimension operator here is (uQu), 
where Q'^'^ is the energy momentum tensor and is the four-velocity of the heat bath and 
it is introduced to restore Lorentz invariance formally in the thermal field theory. In the 
rest frame of the heat bath, we have = (1, 0, 0, 0) which leads to = 1. Note that in our 
calculations, we ignore the heavy quark condensate (ip'ip) since it suppress by inverse powers 
of the heavy quark mass. 

To proceed in calculation of the nonperturbative part, we use the nonperturbative part 
of the quark propagator in an external gluon field, A^(a;) in the Fock-Schwinger gauge, 
x'^A^(x) = 0. In this gauge, the vacuum gluon field A'Jj^{k') is written in terms of gluon field 
strength tensor in momentum space as follows: 

^;:(^') = -^(2vr)^G^,(0)A5(4)(A;'), (18) 

where k' is the gluon momentum. Taking into account one and two gluon lines attached 
to the quark line as shown in Fig. 1, up to terms required for our calculations, the non- 
perturbative part of the temperature-dependent massive quark propagator is obtained as: 

+ (m^ - A{k ■ uf) ^ + A{k- u){k^ - m^) i] (m^B^^m^)}, (19) 
where 0^^ is the traceless gluonic part of the energy- momentum tensor of the QCD. 
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Using the above expression and after straightforward but lengthy calculations, we get the 
following expression for the nonperturbative part: 

/' dx{ ^^^ „ „. „ ^ , . ., [g^(9A(-8g^(-l + x)V-mV(-l + a;)V 

X (-3 - 26a; + 2%x^ - ^x^ + 1x^) + 3m^(5 - 16a; + 26a;2 - 20a;^ + lOa;^) 

+ 2m^g2a;(-6 + 3a; + 16a;^ - 33a;^ + 30a;^ - lOa;^)) + S(-8m^g^(-l + xf 

X a;2(4 + 5a; - \x^ - 2x^ + x^) - q^{-l + xfx^{6 + 19a; - 17a;^ - 4a;^ + 2a;^) 

- 2m^(29 - 92a; + 150a;^ - 116a;=^ + 58a;^) + llm^g^a; 

X (6 - 11a; + Sa;^ + a;^ - 6a;^ + 2x^))) - 2fi(-2g^(-l + x) V 

X (6 - 17x + 19a;2 - Ax^ + 2x^) - m^q\-l + xfx^ 

X (37 - 154a; + 188x^ - Q^x^ + 34a;^) + m^il'd - 64a; + 1023;^ - 76a;^ 

+ 38x^) - 4m^g2a;(-6 + 47x - 116^^ + 143^3 - 102x^ + 34a;^))(g • uf]], (20) 

where A = ^(C^^G"''^) + \{u^QipU^) and B = \{u''QipU^). 



C. Thermal Sum Rules for Physical Quantities 



Now it is time to equate two different representations of the correlation function from 
physical and QCD sides and perform continuum subtraction to suppress the contribution 
of the higher states and continuum. As a result of this procedure we get the following sum 
rule including the temperature-dependent mass and decay constant: 



[m\{T) + Ql] m\{T) 



PS 



-ds 



T-rnp 



(21) 



s2(s + gg) 

where So(^) is temperature-dependent continuum threshold and for simplicity, the temperature- 
dependent width of meson has been neglected. To further suppress the higher states and 
continuum contributions, we also apply the Borel transformation with respect to Qq to both 
sides of the above sum rule. As a result we get. 



fl{T)m\{T)exp{- 



M2- 



ds 



Am? 



M2 



) + BU 



np 
t ) 



(22) 



where the nonperturbative part in Borel scheme is obtained as: 



/■ 

Jo 



dx- 



1 



-exp 



m 



M'^{-l + x)x 



/{3A(8M6(-l + a;)V 



+ m^(l - 2a;)2(l - 2x + 2x^) - m^M^{-\ + x) V 
X (-3 - 2a; + 4^^ - Ax^ + 2x^) + m^M'^xi;^ - lAx + lAx^ 
- 2Ax^ + 8a;^)) - S(3m^(l - 2a;)2(l - 2a; + 2x^) + M\ 
X (6 - 29a; + 31a;^ - 4a;=^ + 2a;^) - rri^M^-l + xfx^(-A - 29x 
+ 43a;^ - 28a;^ + 14a;^) + m^M2a;(8 - 35a; + 22a;^ + 69a;^ - 96a;^ + 32a;^))}. (23) 



13x-^ 
-1 + x)^x^ 
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Here is the Borel mass parameter. Considering Eq. (122|) . the mass squared of the heavy 
axial vector meson alone can be obtained as: 

G{T) 



F{Ty 



where, 



and 



F{T) = I ds 

Jim? 



ptA^) + pas{.s) exp{-—) + 5n7, 



M2- 



d 



(24) 



(25) 



(26) 



D. Numerical Results 



To numerically analyze the sum rules for mass and decay constant, we use the following 
temperature-dependent continuum threshold |l8l |: 



So(T) = So 



+ 4 



T* 



(27) 



where T* = 1.1 = 0.176 GeV with being critical temperature and Sq is the continuum 
threshold at zero temperature. For the temperature-dependent gluon condensate we also 
use 13, B 

(OIG^IO) 



exp 



12(J-1.05 



(28) 



For the thermal average of total energy density (O) we use both results: i) obtained in lattice 
QCD 0,3: 



(6) = 2(05) = 6 X 10-^exp[80(T - 0.1)]{GeV^ 



(29) 



where this parametrization is valid only in the region 0.1 GeV < T < 0.17 GeV. ii) obtained 
via chiral perturbation theory 32l |: 



(9) = (ep + 3 p, (30) 
where (0{^) is trace of the total energy momentum tensor and p is pressure. They are given 



as: 



270 V T 

^/m^T\|/ 15 T 105 \ / mT,\ 
p = 3T(-^) M 1 + + ——^)exp( - -f ) , 



(31) 
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where Ap = 0.275 GeV, = 0.093 GeV and m„ = 0.14 GeV. 

We also use the values = (1.3±0.05) GeV, = (4.7±0.1) GeV and (0 | ^a^G^ \ 0) = 
(0.012 ± 0.004) GeV^ for quarks masses and gluon condensate at zero temperature. Finally, 
we should find the working region for the continuum threshold at zero temperature (sq) 
and Borel mass parameter (M^) such that the physical observables are weakly depend on 
these parameters according to the standard criteria of the QCD sum rules. The continuum 
threshold, Sq is not totally arbitrary and it is correlated to the energy of the first exited 
state of the heavy axial vector meson. Our numerical calculations lead to the intervals 
So = (106 — 110) GeV^ and sq = (15 — 17) GeV^ for the Xbi and xa heavy axial mesons, 
respectively. The working region for the Borel mass parameter is calculated requiring that 
not only the contributions of the higher states and continuum are efficiently suppressed but 
also the contributions of the operators with higher dimensions are ignorable. We get the 
working regions 10 GeV^ < < 35 GeV^ and 5 GeV^ < < 25 GeV^ respectively for 
the Xbi and Xa channels. 
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FIG. 2. Dependence of the mass of the Xbi meson on the Borel parameter at zero temperature. 



Using the above obtained working regions for auxiliary parameters together with the 
other inputs, we plot the dependence on the Borel parameter of the masses and decay 
constants of the heavy axial Xbi and Xci quarkonia at zero temperature in Figs. ([2l|5]). From 
these figures, we see that the results weakly depend on the auxiliary parameters in their 
working regions. The numerical results for the masses and decay constants of the heavy 
axial vector mesons under consideration are depicted in tables I and II. We also compare the 
obtained results with the experimental values in the same tables. From table I we see a good 
consistency of our results with the experimental data. The errors in the results of our work 
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FIG. 3. Dependence of the decay constant of the Xbi meson on the Borel parameter at zero 
temperature. 



4.0- 



3.8- 



3.6 



CD 



3.4- 



3.2- 



3.0- 



10 



10 



15 



20 



25 



s„=17 GeV 
■s„=16 GeV' 
s=15GeV' 



4.0 



3.8 



3.6 



3.4 



3.2 



3.0 



15 

M^{GeV^) 



20 



25 



FIG. 4. Dependence of the mass of the Xci meson on the Borel parameter at zero temperature. 

belong to the uncertainties in calculation of the working regions for auxiliary parameters as 
well as those coming from other inputs. 

At the end of this section we would like to discuss the behavior of the decay constants and 
masses of the heavy axial quarkonia under consideration in terms of temperature. We depict 
the variations of these quantities versus temperature in figures ([Ml]). From these figures, we 
see that the masses and decay constants remain unchanged with the variation of temperature 
up to T = 100 MeV. After this point they start to decrease increasing the temperature. At 
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FIG. 5. Dependence of the decay constant of the Xci meson on the Borel parameter at zero 
temperature. 
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FIG. 6. Dependence of the mass of the Xbi meson on temperature at = 20 GeV^ 



deconfinement or critical temperature, the decay constants decrease about (73-78)%, while 
the masses are decreased about 4%, and 19% for Xbi and Xci states, respectively. The sharp 
decreasing in the values of the decay constants near the deconfinement temperature can be 
considered as a signal for existing the QGP as the new phase of hadronic matter. 
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FIG. 7. Dependence of the decay constant of the Xbi nieson on temperature at = 20 GeV^ 
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FIG. 8. Dependence of the mass of the Xci meson on temperature at = 10 GeV^ 
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m^,, {GeV) 


"^Xbi (GeV) 


Present Work 


3.52 ±0.11 


9.96 ± 0.26 


Experiment [34] 


3.51066 ± 0.00007 


9.89278 ± 0.00026 ± 0.00031 



TABLE I. Values for the masses of the heavy axial Xbi Xci quarkonia at zero temperature. 
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